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CN : ABSTRACT 
O 

\q | Molecular clouds interacting with supernova remnants may be subject to a 

greatly enhanced irradiation by cosmic rays produced at the shocked interface 
between the ejecta and the molecular gas. Over the past decade, broad-band 
observations have provided important clues about these relativistic particles and 
indicate that they may dominate over the locally observed cosmic-ray popula- 
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tion by a significant amount. In this paper, we estimate the enhancement and 
find that the cosmic ray energy density can be up to ~1000 times larger in the 
molecular cloud than in the field. This enhancement can last for a few Myr and 
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leads to a corresponding increase in the ionization fraction, which has important 
consequences for star formation. Ionization fractions in molecular cloud cores 
determine, in part, the rate of ambipolar diffusion, an important process in core 
formation and pre-collapse evolution. Ionization fractions in newly formed cir- 
cumstellar disks affect the magneto-rotational instability mechanism, which in 
turn affects the rate of disk accretion. As estimated here, the increased ioniza- 
tion acts to increase the ambipolar diffusion time by a factor of ~ 30 and thereby 
suppresses star formation. In contrast, the increased ionization fraction reduces 
the sizes of dead zones in accretion disks (by up to an order of magnitude) and 
thus increases disk accretion rates (by a comparable factor). 



Subject headings: Cosmic rays — ISM: Clouds — Magnetohydrodynamics - 
Stars: Formation — Supernovae — Turbulence 
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1. INTRODUCTION 

Although a robust paradigm for star formation within giant molecular cloud complexes 
now exists (Shu et al. 1987), the interplay between stars and the environment from which 
they are born remains poorly understood. This letter considers the interaction between 
supernova remnants (SNRs) and molecular clouds (MCs), and studies its effect on the level of 
cosmic ray fluxes within star forming regions. The cosmic ray flux affects the ionization rate, 
which influences the formation of molecular cloud cores and the evolution of circumstellar 
disks. Specifically, we explore how the cosmic ray production at the SNR/MC boundary 
influences the ionization fraction of the molecular gas. This ionization fraction controls the 
ambipolar diffusion rate (and hence the rate of core formation) and the action of magneto- 
rotational instabilities (MRI) in circumstellar disks (and hence the rate of disk accretion). 

A growing amount of data concerning the interactions of SNRs and MCs has been 
assembled. Huang & Thaddeus (1986) found that about half of SNRs in the outer galaxy 
are spatially coincident with large molecular cloud complexes. More recently, 1720 MHz 
OH maser emission, a signpost for SNR-MC interactions (Wardle & Yusef-Zadeh 2002), has 
been observed from nineteen galactic SNRs - roughly 10% of the sample (Yusef-Zadeh et 
al. 2003). The appearance of this maser emission is important in that it helps constrain the 
densities and temperatures of post-shock gas in these environments. In addition, magnetic 
field strengths and orientations can be deduced from Zeeman splitting and polarization 
studies (Claussen et al. 1997; Brogan & Troland 2001). High resolution X-ray and CO data 
have added to the list of interacting SNR/MC sources (Butt et al. 2001; Dubner et al. 2004; 
Byun et al. 2006; Sasaki et al. 2006). 

Shocks produced by the interaction of SNRs and MCs can transfer ~ 10% of the SN 
explosion energy into cosmic-rays (Dorfi 1999, 2000), and the interaction of these relativistic 
particles with the ambient medium can lead to the production of pions and, subsequently, 
neutral pion-decay gamma-rays (Aharonian et al. 1994). Evidence for this hadronic process 
has been found through the association of at least ten EGRET sources with SNRs expanding 
into MCs (Esposito et al. 1996; Combi et al. 1999, 2001; Torres et al. 2003); five of these 
sources also belong to the aforementioned set of maser SNRs. While the EGRET SNRs have 
gamma-ray luminosities in the 30 MeV - 10 GeV band ranging from L 7 ~ 10 34 to 4 x 10 36 
ergs s -1 , they have not been detected at TeV energies (Rowell et al. 2000; Buckley et al. 
1997). This finding sets important constraints on the distribution of energetic particles in 
these systems and corresponding constraints on particle-acceleration mechanisms in these 
high-density environments (Baring et al. 1999). These SNRs are also typically radio loud 
(F ^ 10 Jy at 1 GHz), perhaps owing in part to synchrotron emission from secondary leptons 
produced via the decay of charged pions (Fatuzzo & Melia 2003). 
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Particle acceleration in the SNR/MC environments can significantly enhance the cosmic- 
ray (CR) density above that of the local background "sea" surrounding these regions (Aha- 
ronian & Atoyan 1996, hereafter AA96; Torres et al. 2003). These energetic particles can 
increase the ionization fractions within the molecular cloud, and thereby influence star for- 
mation. One important process in the formation of molecular cloud cores (stellar birth sites) 
is the removal of magnetic flux through ambipolar diffusion. The diffusion coefficient, and the 
rate, depends on the ionization fraction of the cloud material. As another example, an im- 
portant mechanism for producing accretion in circumstellar disks is the magneto-rotational 
instability (MRI) which depends on the ionization fraction in the disk. The enhanced cosmic 
ray flux considered here produces a corresponding increase in the ionization fraction, which 
in turn affects the rates of these star formation processes. 

2. COSMIC-RAY ENHANCEMENT IN SNR/MC ENVIRONMENTS 

This section estimates the enhancement of cosmic rays in SNR/MC environments (see 
also AA96). The EGRET SNRs have inferred gamma-ray luminosities spanning the range 
L 1 ~ 10 34 to 4 x 10 36 ergs s -1 . This large range in luminosities may be due in part to 
the age distributions of the EGRET SNRs. If produced via a pion-decay mechanism, the 
gamma-ray luminosity initially increases as relativistic particles are continually injected into 
the SNR/MC environment and their population builds up. The luminosity reaches a peak 
once the acceleration mechanism is quenched (Sturner et al. 1996). Indeed, if 10% of the 
SN explosion energy is transferred to relativistic particles (Dorfi 1991, 2000), the gamma-ray 
luminosity of SNR/MC environments would peak at the value 

Ll ~ , M E CX „ 2.4 x 10-W (i) (JL) , (1) 

where a pp = 40 mb is the pp scattering cross-section, rj is the fraction of a relativistic 
proton's energy that goes into the 7To decay photon channel, Eqr is the energy content of the 
relativistic particles, and n p and are, respectively, the proton and molecular hydrogen 
number densities (i.e., n p = 2n# 2 ). 1 The value of 77 can be estimated using the approximation 
for the gamma-ray luminosity due to the decay of neutral pions 

2A°(a) . f dN CR (E~) Jrn 2A°(a) . 

L 7 « {(r pp n p c) / £ 7 — dEj « {a pp n p c) E CR , (2) 
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1 Molecular clouds are highly nonuniform, with densities on different scales ranging over several orders of 
magnitude. However, the mean molecular hydrogen density of MCs is observed to be ~ 10 2 cm~ 2 
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where a is the spectral index of the distribution cINcr/cLE of relativistic protons, and A is 
a dimensionless parameter that depends on a (Crocker et al. 2005). The above equations 
imply i] = 2A°/a, so that r\ = 0.18, 0.10, and 0.60 for a = 2.0, 2.2, and 2.4, respectively. The 
populations of relativistic protons injected into SNR/MC environments are expected to have 
power-law distributions with a = 2 - 2.4 (Jones & Ellison 1991), in good agreement with the 
EGRET observations. Since the diffusion rate depends on the energy of the particles, the 
spectrum of relativistic protons may flatten inside a MC so that Aa « 0.2 (Ballantyne et 
al. 2006). Since the ionization rate depends on the number density of particles (rather than 
the energy density), this change in spectral index implies a corresponding change in number 
density (by a factor of ~ 2), but the exact size of the change depends on the details of the 
diffusion process. In future work, a careful assessment of this effect should be carried out. 

Next we assume that particles build up during an injection phase to their peak values on 
a timescale of ~ 10 4 years, consistent with estimates of the EGRET SNR ages. A fraction of 
these particles will then diffuse into the cloud with which the SNR is interacting (where the 
fraction is determined by the solid angle subtended by the cloud). Diffusion continues until 
the cosmic rays undergo pp scattering with the ambient medium, and thereby experience 
catastrophic losses, on a time scale given by 



where k ~ 0.45 accounts for the inelasticity of pp interactions. Since protons only lose a frac- 
tion of their energy in a given scattering event, the total time of the cosmic ray enhancement 
could be larger than r pp by a factor of 2 or 3. Setting the time scale r pp equal to the diffusion 
time scale r^/ = R 2 /D(E), where D(E) is the energy-dependent diffusion coefficient, one 
can estimate the distance R(E) that particles will diffuse into the cloud before scattering and 
losing their energy. Adopting the expression D(E) = D w (E/10 GeV) 1 / 2 , where D w = 10 28 
cm 2 s^ 1 (AA96), we find 



As noted in AA96, however, values of D 10 ~ 10 26 cm 2 s" 1 may be more applicable in the 
molecular cloud environment (Ormes et al. 1988), leading to an estimate of R(10 GeV) = 
10 pc. This length scale is comparable to the typical sizes of molecular clouds. 

The energy loss rate for a cosmic-ray due to ionization is weakly dependent on energy for 
relativistic protons, and is given by dE/dt ~ 3.6 x 10~ 7 n// 2 eV/s (Mannheim & Schlickeiser 
1994). The corresponding cooling time thus takes the form 




(3) 




10 28 cm 2 s- 1 
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which is considerably longer than the pp collision time for all but the least energetic protons. 
The resulting energy density u cr within the MC during the SNR/MC interaction will be 
u cr ~ 500 eV cm -3 (E C r/10 50 ergs) (i?/10pc)~ 3 . This value is about 10 3 times larger than 
the energy density of the local cosmic-rays (where u cr ~ 0.5 eV cm -3 ). Thus, a significant 
enhancement of cosmic ray energy density can be achieved in molecular clouds through their 
interactions with SNRs. 



3. ENHANCEMENT OF THE IONIZATION FRACTION 

The number density of ionized particles in a MC environment depends on the ionization 
rate (due to cosmic-rays) and a complex recombination process. For steady-state conditions, 
the ionization fraction x = n e /nH 2 is well approximated (Elmegreen 1979; Shu 1983, 1992) 
by the expression 

x ~ 10- 5 (C/lO^s- 1 ) 172 (n/lOW 3 )" 172 . (6) 

This expression is valid for n £ 10 6 cm -3 ((/10~ 17 s -1 ). Assuming that ( scales linearly with 
the cosmic ray energy density, i.e., ( = 10~ 17 s -1 (it cr /0.5 eVcm" 3 ), the ionization fraction 
in molecular clouds interacting with SNRs can be enhanced by up to a factor of ~ 30. This 
scaling is strictly valid when n £ 10 9 cm™ 3 for the enhanced cosmic ray flux considered here. 
The process of ambipolar diffusion occurs primarily at lower densities, and will be fully 
affected by this enhancement. In circumstellar disks, however, the densities are larger near 
the disk midplane, and the full ionization enhancement might not be realized. Notice that 
the ionization rate depends on the number density of cosmic rays, rather than the energy 
density. If the spectrum cINqr/cIE is invariant, the two quantities scale together. However, 
the index a of the spectrum could vary (§2), so that the scaling could differ by a factor of 
~ 4 (which affects the ionization fraction by a factor of 2). 

The timescale required to reach this enhancement level can be estimated by setting 
the ionization rate n e = (n equal to the ratio n e /rj. Using the expression for the ionization 
fraction x given by equation (6), this timescale estimate becomes r 8 = i/( « 3 x 10 4 yr 
(C/lO^s- 1 )- 1 / 2 (n/lO^m- 3 )- 1 / 2 , only about 1000 yr when C ~ 10~ 14 s -1 . The ionization 
enhancement is thus established at virtually the same time as the CR enhancement. Of 
course, this result also means that the ionization enhancement will decay on a comparably 
short timescale, as soon as the CR enhancement is removed due to pp scattering (up to 1 
Myr, as estimated in §2; this time can be extended, as discussed below). 

The discussion so far considers the MC environment to be uniform. However, these MC 
regions are highly nonuniform, exhibiting hierarchical structure that can be characterized in 
terms of dense clumps and cores surrounded by an interclump gas with n ~ 5 — 25 cm" 3 . 
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Clumps have characteristic densities n ~ 10 3 cm -3 and radii r « 0.1 — 1 pc. These dense 
regions occupy a relatively small fraction (2 - 8 %) of the cloud volume, but can account for 
most of its mass (Williams et al. 1995). As a result, cosmic rays propagating through the 
MC environment spend most of their time in a lower density medium (n# 2 ~ 10 cm -3 ), so 
the pp collision timescale would be ~ 3 Myr. Of course, this point is somewhat diminished 
by the fact that these particles also spend a small fraction of their life (~ 2 — 8%) in the 
denser regions. If these particles can survive without undergoing pp scattering for ~ 1 Myr, 
the enhancement in cosmic ray energy density can last for a few Myr. 

The interaction of relativistic protons with the MC medium leads to the production of 
both neutral and charged pions, with each species produced in equal numbers. As a result, 
pp scattering transfers roughly 3% of the initial cosmic-ray energy to secondary electrons 
and positrons - the byproducts of charged pion decays. The injection of these leptons can 
increase the energy density of cosmic-rays above the field by as much as a factor of ~ 30, 
and increases the ionization fraction by a factor of ~ 5 (assuming that ( oc u cr ). The peak 
of the injected secondary lepton distribution occurs at ~ 60 MeV. Near this peak energy, 
leptons traversing neutral atomic hydrogen cool predominantly via electronic excitation on a 
timescale r c ~ 7 x 10 4 (n p /10 2 cm -3 ) -1 yr (Gould 1975). The enhancement in the ionization 
fraction due to secondary leptons would therefore only last for a comparably short time. 



4. EFFECTS ON STAR FORMATION 

The enhancement in cosmic ray flux, which increases ionization levels, can have a sig- 
nificant impact on star formation. Ionization affects the coupling between magnetic fields 
and the (largely neutral) gas. On scales of ~ 0.1 pc, increased ionization leads to greater 
coupling between gas and the magnetic fields, and acts to slow down star formation. On 
smaller scales of ~ 1 AU, increased ionization in circumstellar disks allows for MRI to remain 
active over a greater extent of the disk and thereby allows for larger disk accretion rates. 

1 /2 

In considerations of ambipolar diffusion, the density of ions takes the form p^ = Cpn 
because the volumetric rate of recombinations is proportional to n e rii ~ nf and the volumetric 
rate of ionization is proportional to (n n (Shu 1992). As a result, the constant C oc ^/C- 
Standard considerations of ambipolar diffusion (Mouschovias 1976; Shu 1983, 1992) show 
that the effective diffusion constant is given by D ~ v\j {^Cp) oc C -1 ^ 2 , where va is the Alfven 
speed and 7 is the drag coefficient between ions and neutrals. As the cosmic ray flux increases, 
the diffusion coefficient decreases as the square root of the flux. A 1000-fold increase in the 
cosmic ray flux (§2) thus decreases the effective diffusion constant by a factor of ~ 30, 
and increases the ambipolar diffusion time scale by the same factor. Current observations 
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indicate that the diffusion process is too slow to account for the observed statistics of starless 
molecular cloud cores (e.g., Jijina et al. 1999), so that some mechanism to speed up the 
process is required, even in the absence of supernova enhancement of the cosmic ray flux 
(Ciolek & Basu 2001, Zweibel 2002, Fatuzzo & Adams 2002). The factor of ~ 30 increase 
in timescale due to cosmic ray enhancement essentially shuts down the process. Although 
supernovae have often been invoked as a means to trigger star formation (Elmegreen 1998), 
this mechanism provides a channel for supernovae to inhibit star formation. 

An enhancement in cosmic ray flux also influences disk accretion for star/disk systems 
associated with the interaction region. A growing consensus in the field holds that disk 
accretion is produced by an effective viscosity that is driven by turbulence, which in turn 
is driven by MHD instabilities such as MRI (Balbus & Hawley 1991). In order for MRI to 
operate, and hence for disk accretion to take place, the ionization fraction must be sufficiently 
high so that the gas is well coupled to the field. The inner disk can be ionized by collisions 
(where the number density and temperature are high), and the outer disk can be ionized 
by standard values of the cosmic ray flux, but intermediate regions may have dead zones 
where ionization is too low (Gammie 1996). The enhanced cosmic ray flux acts to increase 
the fraction of the disk that is active, i.e., sufficiently ionized for MRI to operate. 

To illustrate the possible effects of an enhanced cosmic ray flux, consider the disk models 
of Gammie (1996), where the cosmic ray flux has its standard (unenhanced) value. The flux 
is assumed to be ineffective as an ionization source after it falls to 1/e times its initial value. 
The attenuation column density for cosmic ray ionization is S ~ 100 g/cm 2 (Umebayashi & 
Nakano 1981). As a result, only the uppermost 100 g/cm 2 of the disk (and a corresponding 
layer on the other side) experience enough cosmic ray ionization for MRI to operate. With a 
cosmic ray flux that is enhanced by a factor J 7 , the disk will experience cosmic ray ionization 
down to a larger column density £* « (1+ln JF)E . For an optimal enhancement of T « 1000, 
the active column density becomes £* « 8S ~ 800 g/cm 2 (~ 8 times more disk material 
can be ionized). To put this result in perspective, consider a circumstellar disk with mass 
Ma = 0.05M Q , radius = 30 AU, and surface density profile £ oc 1/r. The total surface 
density at 1 AU is thus £ (1 AU) ~ 2400 g/cm 2 . Since the upper and lower 800 g/cm 2 
of the column can be ionized with an enhanced cosmic ray flux, about 2/3 of the disk can 
remain active at r = 1 AU. As a result, the dead zone subtends only the central 33% of 
the disk column density, compared with 92% for no cosmic ray enhancement. Since both 
the disk properties and the cosmic ray flux enhancement will vary from system to system, 
an extensive exploration of parameter space is warranted. Nonetheless, this calculation 
illustrates that supernova enhancements of the cosmic ray flux can have a substantial impact 
on disk accretion. In the most extreme case, the disk accretion rate can be an order of 
magnitude larger than in regions with no cosmic ray enhancement. 



- 8- 



5. CONCLUSIONS 

This letter has shown that a supernova interacting with a molecular cloud can produce 
a substantial enhancement in the flux of cosmic rays within the cloud. These cosmic rays can 
diffuse a distance R ~ 10 pc into the molecular cloud before undergoing pp scattering leading 
to catastrophic energy losses. The resulting flux enhancement is estimated to be a factor of 
~ 10 3 (i?/10 pc) -3 within the diffusion region, and can last for a few Myr after cosmic-ray 
acceleration ceases. The corresponding enhancement in the ionization fraction within the 
molecular cloud is estimated to be ~ 30(i?/10 pc)~ 3 / 2 . This cosmic-ray enhancement may 
have a significant impact on star formation processes, primarily by affecting the ionization 
fraction. The enhanced ionization fraction in MC cores, the sites of protostellar collapse, 
increases the ambipolar diffusion time scale by a factor of ~ 30, so that supernovae inhibit 
further star formation via this mechanism. The enhanced cosmic ray flux also increases the 
ionization levels within circumstellar disks. In this setting, more of the disk can actively 
support MRI, so that supernovae enhancement of cosmic ray flux leads to increased rates of 
disk accretion, by a factor of ~8. 

During the past decade, multiwavelength observations have provided important clues 
about cosmic ray acceleration resulting from the interaction between supernova remnants 
and molecular clouds. While OH maser emission has helped constrain the environment 
of the acceleration sites, EGRET observations have helped constrain the resulting particle 
distributions. Additionally, since the gamma-ray emission above ~ 100 MeV is believed to 
result primarily from the decay of neutral pions produced by the scattering of the cosmic rays 
with the ambient medium, future GLAST observations promise to advance our understanding 
of how cosmic rays diffuse from their acceleration sites and interact with the molecular cloud 
environment. Along with these future observations, additional theoretical work should be 
carried out to explore the details of the diffusion of cosmic rays into molecular clouds, the 
enhanced ionization fractions produced therein, and the subsequent effects on star formation. 
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